Orbital evolution of a test particle around a black hole: higher-order corrections 
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We study the orbital evolution of a radiation-damped binary in the extreme mass ratio limit, 
and the resulting waveforms, to one order beyond what can be obtained using the conservation laws 
approach. The equations of motion are solved perturbatively in the mass ratio (or the corresponding 
parameter in the scalar field toy model), using the self force, for quasi-circular orbits around a 
Schwarzschild black hole. This approach is applied for the scalar model. Higher-order corrections 
yield a phase shift which, if included, may make gravitational-wave astronomy potentially highly 
accurate. 
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The problem of determining the orbital evolution of 
a binary which undergoes radiation damping has be- 
come timely and crucial, as earth-based gravitational- 
wave detectors (LIGO/VIRGO) will become operative 
in the near future, and as a space-based detector (LISA) 
is currently planned to fly in about a decade ||l|. The 
orbital evolution becomes a much more feasible problem 
when M ^ n, M being the mass of the central black 
hole, and /z that of the inspiralling companion. For typi- 
cal parameters of a supermassive black hole at galactical 
centers and a stellar black hole, with M/^ « 10^~^, the 
frequency of the emitted waves is in the good-sensitivity 
band for space-based detectors. LISA is expected to mea- 
sure such waves with amplitude signal-to-noise ratio of 
order 10 — 100 with one year integration time, with 
occurrence rate of a few per year to a few per month . 

For M ^ ^, the companion can be modeled as a point- 
like particle, and the orbital evolution can be studied us- 
ing perturbation theory: fi moves in the fixed background 
spacetime of M; because of the (small) perturbation to 
the metric due to /i, the latter does not move along a 
geodesic of the background. Instead, it moves along an 
accelerated orbit, resulting from its self force (SF). (In an 
alternative equivalent picture, the particle moves along a 
geodesic of a perturbed spacetime.) 

The calculation of the SF for realistic astrophysical sys- 
tems turns out to be rather difficult. Because of that dif- 
ficulty, it is generally hoped that simpler methods, appli- 
cable for simplified systems, may be useful. Specifically, 
an approach dubbed "radiation reaction (RR) without 
radiation reaction forces" (RRWORRF) has been used to 
evolve the orbit and obtain the corresponding waveforms 
[Q, ^. The underlying motivation for that approach is 
roughly the following: So long as the orbital evolution is 
adiabatic, the constants of motion (COM) which describe 
the orbit change only slowly. The rate of change of the 
COM can be inferred from the flux of the correspond- 
ing quantities to infinity and down the event horizon of 
the central black hole, and from this rate of change one 
can infer the orbital evolution. One complication to this 
strategy is that one of the COM, the Carter constant Q, 



is non-additive. It is generally impossible to use balance 
arguments to find Q. For that reason, this approach has 
been used for cases where Q is trivial, such as for circular 
Q or equatorial orbits around a Kerr black hole. The 
conventional approach is to compute the local SF acting 
on fj,, /^^, and infer from it Q. With the knowledge of the 
rates of change of all COM, one would be able to evolve 
the orbit and find the corresponding waveform. That 
approach, however, can only find the orbital evolution 
to leading order in ^/M. [In what follows, we denote 
e — fi/M. For the case in which the particle carries a 
scalar or electric charge q, e = /{y.M).] Consequently, 
conservation laws (and conservation of //) are violated to 
O(e^), and meaningful quantities such as the number of 
cycles A/'cyc spent in a logarithmic interval of frequency /, 
dMcyc/ din f , and the orbit are found to 0(e~^). It can 
be expected that corrections of 0(1) to dMcyc/dlnf and 
the orbit may be of practical importance: when the data 
stream is cross correlated with a theoretical template, 
the cross correlation plummets when the two slip by too 
much. Already when the two lose phase by as little as 
one half-cycle, their overlap integral will be strongly re- 
duced §. High accuracy theoretically-derived templates 
are not necessary for detection of the waves. They are 
needed, nevertheless, for precision gravitational-wave as- 
tronomy. (See also 0.) It is also important to provide 
estimates for the error in the first-order waveforms. 

It is hard to generalize the RRWORRF program to 
handle higher-order corrections. The reason being that 
the RRWORRF program, and indeed any nonlocal ap- 
proach based on conservation laws, ignores conservative 
SF effects, because the fields associated with the latter 
decay faster at infinity than those associated with dissi- 
pation. When integrating over a sphere at infinity, the 
conservative effects are in practice discarded. 

In this paper we propose to use the SF directly to find 
the orbit to 0(1). Dimensional-analysis arguments sug- 
gest that corrections to the orbit-integrated dMcyd d\nf 
are at order unity. We present a detailed analysis of a 
simplified toy model, which allows us to calculate this 
correction accurately. We find that the naive analysis 
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overestimates the effect of interest by a full order of mag- 
nitude. 

For simplicity, we apply the approach for quasi-circular 
equatorial orbits around a Schwarzschild black hole. The 
metric in the usual Schwarzschild coordinates is given by 
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Here, dfl^ — dd^ -t- siv? ■& dcj)^ . Our approach is based 
on finding r{t) and uj(t) perturbatively to O(e^), uj being 
the orbital frequency. Henceforth, we denote by an over- 
dot and a prime (partial) differentiation with respect to 
coordinate time t and r, respectively. 

We use the normalization condition for w", namely 
u°'Ua = —1, to eliminate u* from the equations of mo- 



tion (EOM) Du'/dT = fi- 



,k = t,r,(j), where 



D denotes covariant differentiation compatible with the 
metric (0). We next use the t component of the EOM 
to eliminate w*. We can simplify the EOM to first-order 
(nonlinear) ODEs by taking r = V{r), x — Vx'{r)^ and 
X denotes any quantity. We find the EOM to be 
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where u* = - ZM/r - r^cr - V^/{1 - 2M/r)]. 

Here, a measures the deviation from Kepler's law, i.e., 
= M/r^ + a{r). This relation is gauge independent, 
although each of the terms on the right hand side are 
separately gauge dependent S . We next expand in pow- 



ers of e as a{r) = a 
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, denoting the term in x which is at 
0{e^), and being the self acceleration. We then expand 



the self force as /, 
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Solving perturbatively, we find that 
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Notice, that Eq. (H) is only a formal expression, as we 



do not know //^■'. [Currently, f^''' are known for scalar 
field RR for circular orbits around Schwarzschild . For 
gravitational RR, even fj^^"^ are not known as yet.] This, 
however, is not an important problem: The coefficients 
of the unknown terms in Eq. (|^) are much smaller than 
the coefficients of the other terms, such that their rel- 
ative contribution is small. [From dimensional-analysis 

(2) fll2 

arguments, it can be shown that ~ at{M/ fi)/^ , 

/(^) ^ a,(M/M)/(i)', and /f ~ a^Af/r)^/^- Vi'^', 
and consider their contribution with not much greater 
than unity. In deriving these relations recall that it is 
simplest to analyze the scaling of the four-accelerations, 
and only then to obtain the forces. Notice also that these 
terms vanish at r = 6Af.] We shall henceforth happily 

(2) 

ignore the terms involving . 

We can study the importance of the higher-order 
correction by considering two dimensionless quantities, 
dMcyc/ dilnf) and V/{ruj). Then, we compare these 
quantities between a theoretical template accurate to 
0(e^^), and a template accurate to 0(1). [Modeling the 
actual data stream by the 0(1) template, we can test 
the 0{e~^) template.] V/{ruj) is related to the rate of 
change of the envelope of the chirp wave. Notice that 
rfA/'cyc/ d(ln/) = u)"^ /[{2)iTUj] is gauge independent. The 
difference in these quantities between the 0(1) expres- 
sions and their 0{e~^) counterparts can be expanded us- 
ing the expressions above. We find that 
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Notice that the last two quantities are at 0(1). The 
bracketed factor of 2 should be introduced for scalar field 
RR (and is absent for gravitational RR) . Notice also that 
Eqs. (0) and (||) do not depend on (T(2) . 
The orbit can be obtained by integrating 
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The integrands can be expanded to 0(1) using the ex- 
pressions above. These integrals are easy to do us- 
ing fourth-order Runge-Kutta integration. The triad 
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FIG. 1; The orbit of a scalar charge with q'^/(fiM) = 0.1. The 
upper panel shows the last few orbits before the ISCO. The 
solid (dashed) line is the orbit to 0(1) [0(e~^)]. The lower 
panels display two enlargements of the same orbits: the one on 
the left emphasizes the difference in the r values, and the one 
on the right the difference in phase. Here, x — {r /M) cos (p 
and y = {r/M) smcj). 
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FIG. 2: The magnitude of the higher-order corrections. Upper 
panel: A (^) / {^) as a function of r. The solid line is 
the full effect, the dotted (dashed) line is the contribution of 
the term in Eq. (^ proportional to o"(i) [V(2)]. Lower panel: 
A dMcyc/ d{\n f) as a function of r. The solid line is the full 
effect. The dotted line is the contribution of the term in 
Eq. ^ proportional to and the lower (upper) dashed 

line the contribution of the term proportional to crji) [V(2)]- 



r, t(r), (j)(r) can then be inverted to yield t, r{t), (j){t), from 
which the orbit and the waveform can be reconstructed. 

As currently we do not know fj^^"^ for gravitational RR, 
we specialize next to scalar field RR. [When /-^^^ are ob- 
tained for gravitational RR in any regular gauge, the 
analysis below can be repeated for that case.] Specifi- 
cally, we study the model of a scalar charge q of mass 
fi, which in the absence of self interaction moves along a 
circular and equatorial geodesic around a Schwarzschild 
black hole of mass M. For this model, the local SF was 
computed recently in Refs. The scalar field $ 

satisfies the wave equation □$ = —477/?, where p is the 
scalar charge density. Solving this equation for the ap- 
propriate motion in Schwarzschild spacetime (see Ref. 

for details), we obtain the bare SF by ^'^^^f^j^ = qd^^. 
The bare force t>aroj^F regularized using mode- 

sum regularization, a procedure which yields the physi- 
cal, finite part of the SF, /^^. 

We next choose e — 0.1, and integrate to find the or- 
bit. We started from rinitiai = 30Af and integrated in- 
ward toward the innermost stable circular orbit (ISCO) 
at r = 6M. (We do not consider the correction to the 
location of the ISCO, which is irrelevant to the question 
of interest here.) At intermediate points along the orbit 
we evaluated /^^ and /^^' using a best fit to a smooth 
function. Our results are displayed in Fig. |l| We find 
that the orbit to 0(1) decays slower than the orbit to 
0(e~^). This is indeed expected, as fr^\ which does not 



contribute at 0(e^^), is a repulsive force, which slows 
down the decay of the orbit. Note, that because our cor- 
rection to the orbit is at 0(1), the last few revolutions 
look qualitatively the same also for other choices of e. 

We next study the magnitude of the higher-order ef- 
fect by considering Eqs. (0) and (||). Figure || shows 
A dNcycl rf(ln /) and A (^) / {^f . We first notice that 
both are dominated by the terms in the corresponding 
equations proportional to V(2)- In fact, at the ISCO the 
other terms vanish, such that at the ISCO these quan- 
tities are fully described by the terms proportional to 
V{2)- (Strictly speaking, the ISCO may be defined by 
the requirement that dA/'cyc/ d(ln /) vanishes. As noted 
above, we do not consider here the shift in the location 
of the ISCO, and our previous remarks relate, in fact, to 
r = 6M, not the ISCO proper. Note also that it is im- 
portant to evaluate V(2) ■ the effect is controlled by that 
term at and near r = 6M.) Fitting our results, we find 
that AdAfcyc/ d{lnf) D fi , where D « 5.9. Integrat- 
ing from Tstart — 20M down to r = 6M, we find that 
the 0{e~^) template would slip by about a tenth of a 
cycle compared with the data stream. This magnitude 
is independent of e. V/{rw), however, is at 0(e). We 
thus conclude that changes in amplitude will generally 
be unimportant for small values of e. 

Finally, we present in Figure || the waveforms to 0(1) 
and to 0(e~^). We estimate the waveforms by means of 
the usual "restricted waveform" approximation ||ll|, i.e., 

h{t) = C vl{t)cos[(j) 

wave it)]: (10) 
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FIG. 3: The waveforms. In all panels the solid (dashed) line 
describes the waveform to 0(1) [0(e~^)]. The upper left panel 
describes the waves from an early part of the orbit. The lower 
panel describes the last few cycles before the ISCO. The upper 
right panel is a magnification of the last maximum. 

where v^^ — {dcj)/ dty^'^, and C is an amplitude coeffi- 
cient which depends on the distance of the source of the 
waves. (For gravitational waves 0wavc = 20, and that 
for scalar field waves (jiwave = 0-) Notice that our results 
above are independent of our adoption of the restricted 
waveform. The two waveforms are initially in phase. A 
phase difference, described by Eq. (^ builds up slowly, 
and is clearly visible in the last few cycles. In addition 
a difference in amplitude is also clearly seen. The latter 
effect, however, is an artifact of the relatively large value 
of e. This change in amplitude decreases linearly with e. 
The phase shift, on the other hand, is independent of e. 

We emphasize that our conclusions are specific to 
scalar field RR. When gravitational RR is considered, the 
magnitude of the higher-order effects could be different. 
For gravitational RR we pick up an obvious factor of 2 in 
Eq. ^ , such that one may expect the higher-order effect 
to be even more important for gravitational RR. Deter- 
mination of ^ (and its gradients) for gravitational RR 
is required in order to quantify this effect accurately. 

An important issue is the possibility to use templates 
to 0{e~^), to cross correlate against the data stream. 
One may hope, that templates to 0(e^^) for slightly dif- 
ferent parameters may be a good approximation, such 
that one does not have to build templates to 0(1). We 
believe that such an approach may indeed be useful for 
search purposes. This is in general impossible to do to 
high accuracy when the cross correlation is done with 
the full wave train: The waveform is a power series in 
e, in which there are different functions of r at each or- 
der . ^Wh£n_one_clmnges_£_to_^^ 



Indeed, numerical experiments suggest to us that when 
the parameter space is one dimensional, such techniques 
are not very useful. [In the scalar field toy model we 
have, in fact, a two-dimensional parameter space, which 
includes and n/M, which can be varied indepen- 
dently. In the more realistic case of gravitational RR, 
however, q/ n, where q is loosely a charge for the grav- 
itational field of fi (or an active mass), is determined 
uniquely as unity by the Equivalence Principle. Moti- 
vated by the gravitational RR problem, we considered 
an effective one-dimensional parameter space also for the 
scalar field toy model, and varied only q^ /{fiM).] On the 
other hand, when more complicated orbits are concerned 
(e.g., non-equatorial orbits, when both the particle and 
the central black hole are spinning), the parameter space 
will be significantly larger, such that such an approach 
may be quite successful. Because of the very slow in- 
crease in the phase difference between the 0(e~^) tem- 
plate and the data stream, it appears to us that data 
analysis techniques based on dividing the wave train into 
many chunks, in which the change in phasing is negligi- 
ble, may be very successful for this problem. However, 
we emphasize that for realization of the potentially high 
accuracy gravitational-wave astronomy, one would have 
to take higher-order corrections, such as those considered 
here, into consideration. 
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